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Harrison's tight-binding tlieory provides an excellent qualitative description of the electronic struc- 
ture of the elements across the periodic table. However, the resulting band structures are in signifi- 
cant disagreement with those found by standard methods, particularly for the transition metals. For 
these systems we developed a new procedure to generate both the prefactors of Harrison's hopping 
parameters and the onsite energies. Our approach gives an impressive improvement and puts Harri- 
son's theory on a quantitative basis. Our method retains the most attractive aspect of the theory, in 
using a revised set of universal prefactors for the hopping integrals. In addition, a new form of onsite 
parameters allows us to describe the lattice constant dependence of the bands and the total energy, 
predicting the correct ground state for all transition, alkaline earth and noble metals. This work 
represents not only a useful computational tool but also an important pedagogical enhancement for 
Harrison's books. 



I. INTRODUCTION 

Walter Harrison developed an elegant analytic 
theory of the electronic structure of solidsiS . This 
theory has been very sucessful in providing a physi- 
cal understanding of the electronic structure and the 
characteristics of bonding. However, Harrison's the- 
ory of solid state has limited ability to produce accu- 
rate numerical results for the band structure, density 
of states and the relative stability of different crystal 
structures. 

In this work, we have set out to put Harrison's ap- 
proach on a quantitative foundation. We have now 
realized that it is possible to put the Slater-Koster 
parameters in the form given by Harrison but with 
new prefactors and determine new onsite parameters. 
The result is that we retain the universality of Har- 
rison's parameters, which means the same prefactors 
for all transition, alkaline earth and noble metals, but 
with different onsite terms for each element. It is 
clear to us that this approach, perhaps slightly modi- 
fied, may be extended to cover the rest of the periodic 
table. We have also succeeded with a small number 
of additional parameters to describe the volume and 
structure dependence of the energy bands and, there- 
fore, obtain total energies and predictions of relative 
stability. 

Harrison has opted for simplicity in the LCAO 
approach and has created a set of universal hopping 
parameters that can easily be used to perform cal- 
culations. In the tables of his books, Harrison uses 
atomic energies as onsite parameters in his Hamilto- 
nians, which is the main shortcoming of directly using 
these tables, to perform sufficiently accurate band 
structure calculations. However, Harrison pointed 
out(^,p561) that atomic values for the d-state ener- 
gies need to be corrected for differences in d-state 
occupancy, and gave a way for doing that in the case 
of Cr. 



ues of the onsite parameters for the transition metals 
Nb and Pd by using Harrison's hopping parameters 
and uncorrected atomic term values. We compared 
the results of a 6 x 6 Harrison Hamiltonian(without 
p orbitals) as given in Harrison's book and we found 
that the energy bands created this way are in serious 
disagreement with Augmented Plane Wave(APW) 
results (see left Fig. 01. We also tested a 9 x 9 (with 
p orbitals) Harrison Hamiltonian with all hopping 
prefactors kept at Harrison's values, but the onsite 
parameters modified by fitting the energy bands to 
APW calculations.^'"* This modification gave us bet- 
ter results in the d— bands, but there was still a large 
error for the s-like first band(see right Fig. 0). Our 
conclusion is that Harrison's theory can only give a 
qualitative description of the band structure of the 
transition metals even if we fit the onsite terms to 
first-principles results. 
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FIG. 1: Harrison's energy bands of Nb compared to 
APW. Left graph created by 6 x 6 Hamiltonian; right 
graph created by 9 x 9 Hamiltonian and fitted onsite en- 
ergies. 



We illustrate the importance of correcting the val- 
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II. ENERGY BANDS AND DENSITY OF 
STATES 

We have developed a procedure that while main- 
taining the simplicity of Harrison's approach gives 
an impressive improvement that puts the theory on 
a quantitative basis. To accomplish this we have 
made the following modifications to Harrison's the- 
ory: (1) We introduced a p onsite energy as an addi- 
tional parameter to the s and d onsite energies used 
by Harrison, and fit them all to APW results. (2) 
We modified the sp hopping integrals of Harrison, by 
introducing a dimensionless parameter 7^ as follows: 



Vu 



mi' 



(1) 



The parameter 7s provides more flexibility to fit the 
first and sixth bands. (3) We obtained new hopping 
prefactors by simultaneously fitting the APW energy 
bands of the following 12 transition metals: V, Cr, 
Ni, Cu, Nb, Mo, Pd, Ag, Ta, W, Pt and Au. In this 
fit, all 12 elements have the same common prefactors 
"Hii'm but different(for each element) onsite energies s, 
p and d, and also different values for the parameters 
7s and that appear in the hopping parameters. 
Our Hamiltonian corresponds to an orthognal basis 
set as in Harrison. We did the above fitting at the 
equilibrium lattice constants of the structure, which 
is the ground state of each element, and included 
interactions of nearest, second-nearest, and third- 
nearest neighbors for the hcc structure and nearest 
and second-nearest for the fee structure. Of course, 
using more neighbors than Harrison did in the fit 
would automatically make changes in the parame- 
ters, even if one were fitting the same bands. Us- 
ing the parameters determined with the above proce- 
dure, we reproduced APW energy bands and density 
of states(DOS) remarkably well, not only for the 12 
elements originally fitted, but also for the rest of the 
transition metals, the alkaline earth and the noble 
metals, as seen in Fig. |2 and Fig. O for four of the 
elements. Our new Hamiltonian prefactors, common 
for all metals, together with Harrison's original pref- 
actors are shown in Table The onsite terms and 
the parameters 7^ and for each element are shown 
in Table ini 



TABLE I: Harrison's hopping prefactors and our modi- 
fied values. 
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FIG. 2: APW and Modified Harrison's energy bands of 
Ag, Pd, Nb and Mo. The solid line is the modified Har- 
rison result and dash-dotted line is the APW result. 
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FIG. 3: APW and Modified Harrison's density of states 
of Ag, Pd, Nb and Mo. 
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TABLE II: Onsite parameters, 7s and rd- 
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"Ferromagnetic spin up 
'Ferromagnetic spin down 



TABLE III: Magnetic Moments of Fe, Co and Ni. 
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FIG. 4: Energy bands of hep Ti and Ru. 



III. TOTAL ENERGY 

Next we address the issue of fitting total en- 
ergy results. In order to do this we follow the 
Naval Research Laboratory tight-binding(NRL-TB) 
methodology"^, which uses parameters which are 
transferable between structures^. In most TB ap- 
proaches as well as in all the so-called "glue" poten- 
tial atomistic methods one writes the total energy as 
a sum of a band energy term( sum of eigenvalues) 
and a repulsive potential G[n{r)] that can be viewed 
as replacing all the charge density n{r) dependent 
terms appearing in the total energy expression of the 
density functional theory. The NRL-TB method has 
the unique feature that eliminates G by the following 
ansatz. We define a quantity Vq: 



We used the prefactors of Table |l] with new on- 
site energies 7s and to fit the rest of the transition 
metals, including those with hep ground states. For 
the hep metals, we fitted energy bands of fee struc- 
tures at the equilibrium lattice, and found that our 
parameters produce good transferability, ie. repro- 
duced the hep energy bands very well without fitting 
them. The hep energy bands of Ti and Ru are shown 
in Fig. ^ We also fitted energy bands of the fer- 
romagnetic elements Fe, Co and Ni, and calculated 
magnetic moments of the three elements at the ex- 
perimental lattice constant. Table IIIII shows good 
agreement of magnetic moments of Fe, Co and Ni 
with experimental values. 
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(2) 



where is the number of valence electrons. 
We then shift all first-principles eigenvalues ei{k) 
by the constant Vq and define a shifted eigenvalue: 



e'^{k) = e,{k) + Vo 



(3) 



The results of this manipulation is that the first- 
principles total energy E is given by the expression: 



(4) 



We note the constant Vo is different for each vol- 
ume and structure of the first-principles database. 
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The reader should recognize that we have shifted 
each band structure by a constant, retaining the ex- 
act shape of the first-principles bands. It should also 
be stressed that all this is done to the first-principles 
database before we proceed with the fit that will gen- 
erate the TB Hamiltonian. In our trearment of fer- 
romagnetic systems the total energy is equal to the 
sum of spin up and spin down shifted eigenvalues. 
The difference of these two sums could be viewed as 
representing the exchange energy. 

We write the onsite energies in a polynomial form: 



hiiipi) = a/ + Pipl + "fipl + Sipf 



(5) 



where I is an angular momentum index, and pi is 
an atomic-like density that has the form: 



p,=Y,e^P[->^^R^3]Fc{R^J) 



(6) 



where, i?^ is the distance between atom i and j , 
and Fc is a smooth cut-off function that was used to 
limit the range of parameters^ 



Fm = (1 



(7) 



We take i?o to be in the range of lO.Oao ~ 14.0ao, 
and Ri — 0.5ao ( where ao is Bohr radius), which 
effectively zeros all interactions for neighbors more 
than 14.0ao apart. Typically, depending on the struc- 
ture and lattice constant, this cut-off function will 
include 50 ~ 80 neighboring atoms. 

The parameters A, a;, 7; and (5; are determined 
by fitting total energies following the NRL-TB proce- 
dure as stated above. The hopping parameters were 
calculated using the modified prefactors of Table HJ 

We fitted total energies of all transition metals 
to the APW results^ at several lattice constants of 
bcc, fee and sc structures. We successfully repro- 
duced the ground-state, the order of crystal struc- 
tures and the bulk modulus. Our parameters also 
place the energies of hep and sc structures, which we 
did not fit, at reasonable values. As an example, we 
present energy-volume relationships for four transi- 
tion metals in Fig. |S1 and for the hep metals Ti and 
Co in Fig. El which again show the correct ordering 
of crystal structures. We also present in Table IIVI 
the equilibrium lattice constants and bulk moduli of 
transition metals. 

An inspection of Table IIVI reveals that our ap- 
proach matches very well the LDA lattice constants 
underestimating the experimental values by 1 — 2% 
for the transition metals and by 4—5% for the alkaline 
earths. The bulk moduli have larger deviation from 
experiment as is usually the case in the LDA. For 
the hep metals both the lattice parameters and bulk 
moduli are also within the LDA predictions except 
for Tc, Os and Y. Those results could be improved if 
we include the hep lattice in the fitting database. 
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FIG. 5: Total energies of Ag, Pd, Nb and Mo. The 
solid line denotes bcc APW results and dash-dotted line 
denotes fee APW results. The diamond, square, triangle 
and X symbols represent bcc, fee, hep and sc TB results 
respectively. 
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FIG. 6: TB total energies of Ti and Co. The solid line, 
dotted line and dashed line represent hep, bcc and fee 
total energies respectively. The symbol -I- denotes simple 
cubic total energy. In the right graph, the thick lines 
and thin lines represent ferromagnetic and paramagnetic 
Co, and the symbols X and -I- show ferromagnetic and 
paramagnetic simple cubic respectively. 



IV. CONCLUSION 

To recapitulate, we have accomplished two goals. 
In the first we have reevaluated the ten universal pref- 
actors in Harrison's hopping parameters and redeter- 
mined the s,p,d onsite energies together with the pa- 
rameters 7s and r^. This enables us to calculate very 
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TABLE IV: Equilibrium lattice eonstants and bulk mod- 
uli for the experimentally^ observed ground-state struc- 
tures of the elements, comparing the results of first- 
principles and tight-binding parametrization results. 
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accurately the band structure of all the transition, 
alkaline earth and noble metals. For the second goal 
we have used a polynomial form for the onsite ener- 
gies which, with the addition of 15 new parameters, 
provides a total energy capability for our Modified 
Harrison theory. 



Finally, we wish to stress that this work consti- 
tutes not only an efficient computational method but 
also a valuable addendum to Harrison's books. 
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